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Abstract 

A resonance search has been made in the K^p and K^p invariant-mass spectrum 
measured with the ZEUS detector at HERA using an integrated luminosity of 121 
pb~^. The search was performed in the central rapidity region of inclusive deep 
inelastic scattering at an ep centre-of-mass energy of 300-318 GeV for exchanged 
photon virtuality, Q^, above 1 GeV^. Recent results from fixed-target experi- 
ments give evidence for a narrow baryon resonance decaying to K^n and p, 
interpreted as a pentaquark. The results presented here support the existence 
of such state, with a mass of 1521.5 ± 1.5(stat.);*^f ^(syst.) MeV and a Gaussian 
width consistent with the experimental resolution of 2 MeV. The signal is visible 
at high and, for > 20 GeV^, contains 221 ± 48 events. The probability of 
a similar signal anywhere in the range 1500-1560 MeV arising from fluctuations 
of the background is below 6 x 10~^. 
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1 Introduction 



Recent results from fixed-target experiments give evidence for the existence of a narrow 
baryon resonance with a mass of approximately 1530 MeV and positive strangeness [1], 
seen in the K^n decay channel. These results have triggered new interest in baryon 
spectroscopy since this baryon is manifestly exotic; it cannot be composed of three quarks, 
but may be explained as a bound state of five quarks, i.e. as a pentaquark, = uudds. 
A narrow baryonic resonance close to the observed mass is predicted in the chiral soliton 
model [2]. The quantum numbers of this state also permit decays to p and p 
(denoted as K^g p (p))- Evidence for a corresponding signal has been seen [3] in this 
channel by other experiments. Evidence for two other pentaquark states has also been 
reported recently [4,5]. 

The Particle Data Group (PDG) [6] lists a number of 'S bumps', unestablished resonances 
observed with low significance by previous fixed-target experiments. The possible presence 
of these resonances in the mass region close to the production threshold of the Kg p (p) 
final state complicates the search for pentaquarks in this decay channel. 

The 0^ state and the S bumps discussed above have never been observed in high- 
energy experiments, where hadron production is dominated by fragmentation. This paper 
presents the results of a search for narrow states in the Kgp[p) decay channel in the cen- 
tral rapidity region of high-energy ep collisions, where particle production is not expected 
to be influenced by the baryon number in the initial state. The analysis was performed 
using deep inelastic scattering (DIS) events measured with exchanged-photon virtuality 
> 1 GeV^. The data sample, collected with the ZEUS detector at HERA, corresponds 
to an integrated luminosity of 121 pb~^, taken between 1996 and 2000. This sample is 
the sum of 38 pb~^ of e~^p data taken at a centre-of-mass energy of 300 GeV and 68 pb"^ 
taken at 318 GeV, plus 16 pb~^ of e~p data taken at 318 GeV. 



2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [7]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [8], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
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polar-angle^ region 15° < 6* < 164°. The transverse-momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pr © 0.0065 © 0.0014/pr, with pr in GeV. To estimate the 
ionization energy loss per unit length, dE/dx, of particles in the CTD [9], the truncated 
mean of the anode-wire pulse heights was calculated, which removes the lowest 10% and 
at least the highest 30% depending on the number of saturated hits. The measured dE/dx 
values were corrected by normalising to the average dE/dx for tracks around the region of 
minimum ionisation for pions, 0.3 < p < 0.4 GeV. Henceforth, (iE'/c/x is quoted in units 
of minimum ionising particles (mips). The dE/dx distribution for electrons has a roughly 
Gaussian shape centred about dE/dx ~ 1.4 mips with width 0.14 mips, corresponding 
to a resolution of ~ 10%. 

The high-resolution uranium-scintillator calorimeter (GAL) [10] consists of three parts: 
the forward, the barrel and the rear calorimeters. The smallest subdivision of the calorime- 
ter is called a cell. The GAL energy resolutions, as measured under test-beam conditions, 
are a{E)/E = 0.18/ ^/E for electrons and a{E)/E = 0.35/V^ for hadrons, with E in 
GeV. A presampler [11] mounted in front of the calorimeter was used to correct the energy 
of the scattered electron^. The position of electrons scattered close to the electron beam 
direction is determined by a scintillator-strip detector [12]. 

The luminosity was measured using the bremsstrahlung process ep epy with the 
luminosity monitor [13], a lead-scintillator calorimeter placed in the HERA tunnel at 
Z = -107 m. 

3 Event simulation 

Inclusive DIS events were generated using the Ariadne 4.08 Monte Garlo (MG) model [14] 
interfaced with Heracles 4.5.2 [15] via the Djangoh 1.1 program [16] in order to 
incorporate first-order electroweak corrections. The Ariadne program uses the Lund 
string model [17] for hadronisation, as implemented in Jetset 7.4 [18]. 

Before detector simulation, the Kg p (p) invariant-mass distribution was calculated from 
the true Kg and (anti)protons in the mass range up to 1700 MeV. No peaks were found, 
indicating that no reflection from known decays are expected to generate a narrow peak 
in these decay channels. 

The generated events were passed through a full simulation of the detector using Geant 
3.13 [19] and processed with the same reconstruction program as used for the data. The 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 

^ Henceforth the term electron is used to refer both to electrons and positrons. 
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detector-level MC samples were then selected in the same way as the data. The generated 
MC statistics were about three times higher than those of the data. 



4 Event sample 

The search was performed using DIS events with > I GeV^, the largest event sample 
for which no explicit trigger requirement was imposed on the hadronic final state. A 
three-level trigger [7] was used to select events online. At the third level, an electron with 
an energy greater than 4GeV and a position outside a box of 24 x 12 cm^ on the face 
of the calorimeter was required. The trigger has a high acceptance for ^ 2 GeV^. 
However, data below ^ 20 GeV^ are strongly affected by prescales which were applied 
to the inclusive triggers to control data rates. 

The Bjorken scaling variables x and y, as well as Q^, were reconstructed using the electron 
method (denoted by the subscript e), which uses measurements of the energy and angle 
of the scattered electron, or using the Jacquet-Blondel (JB) method [20]. The scattered- 
electron candidate was identified from the pattern of energy deposits in the GAL [21]. 

The following requirements were imposed: 

• Ql>l GeV^; 

• -Eg' > 8.5 GeV, where E^^ is the corrected energy of the scattered electron measured 
in the GAL; 

• 35 < 5 < 60 GeV, where 6 = X^-Ej(l — cos^j), Ei is the energy of the ith calorimeter 
cell, 6i is its polar angle and the sum runs over all cells; 

• ye< 0.95 and yjB > 0.01; 

• I ^vertex |< 50 cui, whcre ^vertex is the vertex position determined from the tracks. 

The present analysis is based on charged tracks measured in the GTD. The tracks were 
required to pass through at least five GTD superlayers and to have transverse momenta 
Pt ^ 0.15 GeV and pseudorapidity in the laboratory frame |?7| < 1.75, restricting the 
study to a region where the GTD track acceptance and resolution are high. Gandidates 
for long-lived neutral strange hadron decaying to two charged particles are identified by 
selecting pairs of oppositely charged tracks, fitted to a displaced secondary vertex. Events 
were required to have at least one such candidate. 

After these selection cuts, a sample of 1 600 000 events remained. 
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5 Reconstruction of Kq candidates 



The Kg mesons were identified by their charged-decay mode, Kg tt+tt^. Both tracks 
were assigned the mass of the charged pion and the invariant mass, M{'k~^tt~), of each 
track pair was calculated. The Kg candidates were selected by imposing the following 
requirements: 

• M(e"'"e~) > 50MeV, where the electron mass was assigned to each track, to ehminate 
tracks from photon conversions; 

• M(p7r) > 1121 MeV, where the proton mass was assigned to the track with higher 
momentum, to eliminate A and A contamination of the Kg signal; 

• 483 < M(7r+7r-) < 513 MeV; 

• Pt{K'^s) > 0.3 GeV and < 1.5. 

Figure H shows the invariant-mass distribution for Kg candidates for > 1 GeV^ . A fit 
using two Gaussian functions plus a first-order polynomial function was used. The number 
of Kg candidates was 866800 ± 1000, with a background under the peak constituting 
approximately 6% of the total number of candidates. The peak position was m^^o = 
498.12 ±0.01(stat.) MeV, which agrees with the PDG value of 497.67±0.03 [6] within the 
calibration uncertainty of the CTD. The width is consistent with the detector resolution. 

6 Selection of p{p) candidates 

The (anti)proton candidate selection used the energy-loss measurement in the CTD, 
dE/dx. Figure 121 shows the dE/dx distribution as a function of the track momentum 
for positive and negative tracks. Tracks fitted to the primary vertex were used, with the 
exception of the scattered-electron track. Tracks were selected as described in Section |3] 
In addition, only tracks with more than 40 CTD hits were used to ensure a good dE/dx 
measurement. The tracks were then selected by requiring / < dE/dx < F, where / 
and F, motivated by the Bethe-Bloch equation, are the functions: / = 0.35/p^ + 0.8, 
F = 1.0/p2 + 1.2 (for positive tracks) and / = 0.3// + 0.8, F = 0.75/p2 + 1.2 (for neg- 
ative tracks), where p is the total track momentum in GeV. These cuts were found from 
an examination of dE/dx as a function of p and were checked by studying (anti)proton 
candidate tracks from A (A) decays. The proton band was found to be broader than that 
of the antiproton. There is also a clear deuteron band for positive tracks, which suggests 
a small contribution from secondary interactions. To remove the region where the proton 
band completely overlaps the pion band, the proton momentum was required to be less 
than 1.5 GeV. Finally, a cut requiring dE/dx > 1.15 mips was applied. 
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After these cuts, the purity of the proton sample, estimated from the MC simulation, 
is around 60%, varying from 96% at low momentum to 17% at the highest accepted 
momenta. Applying a higher dE/dx cut leads to higher purity, but reduces the acceptance 
for protons in the high-momentum region and reduces the statistics. 

7 Reconstruction of Kg p (p) invariant mass 

The Kg p (p) invariant mass was obtained by combining Kg and (anti)proton candidates 
selected as described above, and fixing the Kg mass to the PDG value. The resolution 
of the Kg p (p) invariant-mass measurement was estimated using MC simulations to be 
2.0 ± 0.5 MeV in the region near 1530 MeV, for both the Kg p and the Kg p channels. 

The resolution was independently verified from the data by reconstructing the K*{892) 
peak, assuming that the momentum and angle resolution for (anti)protons is similar to 
that for pions. Tracks passing the same angle and momentum cuts as the proton candidate 
sample were taken from the pion band, assigned the pion mass, and combined with the 
Kg candidates. A fit was performed using a Breit-Wigner function convoluted with a 
Gaussian to describe the resolution. The width of the Breit-Wigner function was fixed 
to the natural width of the K* resonance of 50.8 MeV [6]. The resulting width of the 
Gaussian was 5 ± 2(stat.) MeV, independent of the pion charge. 

8 Results 

The Kgp (p) mass spectrum, in the range from threshold to 1700 MeV, is shown in Fig.|21 
for various regions of the DIS phase space. Figures O^i'-Eli show the mass distribution 
integrated above a minimum Q ranging from 1 GeV^ to 50 GeVl The data show signs 
of structure below about 1600 MeV. For > lOGeV^, a peak is seen in the mass 
distribution around 1520 MeV. Figures andEF show the > 1 GeV^ sample divided 
into two bins of the photon-proton centre-of-mass energy, W. A peak is seen in the lower 
W bin. 

The expectation from the MC simulation, scaled to agree with the data in the mass region 
above 1650 MeV, is also shown. If normalised to the luminosity, the simulation lies below 
the data by a factor of approximately two (not shown). Even after scaling, the data are 
not well described by the simulation for masses below 1600 MeV, and no structure is seen 
in the simulated data. The PDG reports a S bump at 1480 MeV, and several states above 
1550 MeV. None of these states are included in the simulation, which includes only well 
established resonances. 
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Several functional forms were fit to the data for > 20 GeV^ over the mass range from 
threshold up to 1700 MeV to estimate the significance of the peak, as well as its position 
and width. A background of the form 

Pi(M -rrip- rriKo)''^ x (1 + P^iM - - m^o)), 

where M is the Kgp (p) candidate mass, rrip and m^o are the masses of the proton and the 
Kg, respectively, and Pi, P2 and P3 are parameters, was found to give a good description 
of the data when combined with Gaussians to describe the signal near 1520 MeV as 
well as possible contributions from S bumps. The parameters of the Gaussians and the 
background were left free. A bin-by-bin minimisation was used. The results of the fit 
using the background function plus one and two Gaussians are shown in Tabled Reducing 
the number of parameters in the background function significantly reduces the quality of 
the fit. Adding a third Gaussian does not significantly improve the fit. The result of the 
fit using two Gaussians is shown in Fig. |3] The second Gaussian significantly improves 
the fit in the low mass region, and has a mass of 1465.1 ± 2.9(stat.) MeV and a width 
of 15.5 ± 3.4(stat.) MeV and may correspond to the S(1480). However, the parameters 
and significance of any state in this region are difficult to estimate due to the steeply 
falling background close to threshold. The signal peak position is 1521.5 ± 1.5(stat.) MeV, 
with a measured Gaussian width 6.1 ± 1.6(stat.) MeV, consistent with the resolution. The 
fit gives 221 ± 48 events above the background, corresponding to 4. 60". The equivalent 
estimate for the single Gaussian fit is 3.9a. The P^r degree of freedom for the fit 
over the whole fitted range of masses is 35/44. Over the same range, the P^r degree 
of freedom for the fit with no Gaussians is 69/50, which is an acceptable fit. However, 
this value is dominated by contributions from the high-mass region. A number of MC 
experiments were carried out, using the background fit with zero gaussians as a starting 
distribution and generating random data using Poisson statistics. The probability of a 
fiuctuation leading to a signal with 3.9cr significance or more in the mass range 1500-1560 
MeV and with a Gaussian width in the range 1.5-12 MeV, was found to be 6 x 10"^. The 
same exercise was carried out using the background plus the 1465 MeV Gaussian as the 
starting distribution, and the probability was found to be about a factor of ten lower. 

The Gaussian function was replaced by a Breit-Wigner function convoluted with a Gaus- 
sian to describe the peak near 1522 MeV. The width of the Gaussian distribution was 
fixed to the experimental resolution to obtain an estimate of the intrinsic width of the 
signal. The extracted width was F = 8 ± 4(stat.) MeV. 

The fit was repeated for different values of the minimum cut. Above ~ 10 GeV^ 
both the signal and background are consistent with having a dependence, similar 

to the inclusive cross section. A MC study in which E"*" particles were modified to have a 
mass of 1530 MeV and artificially forced to decay to Kgp{p) indicated that, at low Q^, the 
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impact of detector acceptance on the visible cross section is important; for < 10 GeV^ 
the number of selected candidates rises more slowly than as well as more slowly 

than the background. Such acceptance effects may be the reason for the absence of a 
clear signal at low and high W. However, this suppression of the signal may also be 
related to the unknown production mechanism of the signal. 

The invariant-mass spectrum was investigated for the Kg p and Kg p samples separately. 
The result is shown as an inset in Fig.|3]for > 20 GeV^, compared to the fit to the com- 
bined sample scaled by a factor of 0.5. The results for two decay channels are compatible, 
though the number of Kg p candidates is systematically lower. The mass distributions 
were fitted using the same function as the combined sample and gave statistically con- 
sistent results for the peak position and width (not shown). The number of events in 
the Kgp channel is 96 ± 34. If the signal corresponds to the G"^, this provides the first 
evidence for its antiparticle. 

If an isotensor state is responsible for the signal, a B"*""*" signal might be expected in the 
K~^p spectrum [22]. The K^^p^K^p) invariant mass spectra were investigated for a wide 
range of minimum values, identifying proton and charged kaon candidates using dE/dx 
in a kinematic region similar to that used in the K^ p [p) analysis. No peak was observed 
in the K^p spectrum, while a clean lOa signal^ was observed in the K~p spectrum 
at 1518.5 ± 0.6(stat.) MeV, corresponding to the A(1520)Do3- Performing a fit using a 
Gaussian fixed to the detector resolution convoluted with a Breit-Wigner gives an intrinsic 
width of 13.7 ± 2.1(stat.) MeV, consistent with the PDG value of 15.6 ± 1.0 MeV [6]. 

8.1 Systematic checks 

A number of checks have been carried out to study possible reflections from known states 
and to verify the robustness of the 1522 MeV peak. 

Tracks from the proton band within twice the width of the K* peak were removed. Ac- 
cording to a MC calculation, this cut increased the purity of the proton sample by 15% 
and reduced the statistics by a factor of 2.5. The resulting peak position was unchanged. 

The energy of the proton candidates was required to be higher than that of the i^^, to 
reduce the combinatorial background [23]. Using this cut and the i^*-rejection cut, a 
peak may be seen in the mass spectrum even in the > \ GeV^ sample. However, this 
combination of cuts leads to a complicated background shape, making the significance 
and the mass of the signal difficult to evaluate. 

^ The signal for A(1520) in the K~^p spectrum has the same number of events and significance as the 
A(1520) signal in the K^p spectrum. 
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The invariant mass was reconstructed from the data using the mass hypoth- 

esis for proton candidates, to see whether Ds decays contribute to the signal. It was 
verified that heavy- flavour particles cannot contribute to this spectrum if M[K^p{p)) < 
2000 MeV. 

The Kg candidates were combined with primary tracks in the region dE/dx < 1.15 mips 
and p < 1.5 GeV, where pions are expected to dominate over (anti)protons. The invariant 
mass was reconstructed using the same procedure as before, applying the proton-mass 
hypothesis for selected tracks. In addition, the mass distribution was calculated using 
the pion band of the dE/dx plot to select proton candidates, and by using proton and 
Kg candidates from different events. In none of these cases was any structure seen in the 
mass distribution. 

The robustness of the peak was also checked by varying the event and track selections. 
The maximum momentum for the proton candidates was changed in a wide range from 
1.2 GeV to 4.0 GeV. The dE/dx cut was varied in the range 1.1 to 1.3 mips. No change 
was seen in the peak position. 

8.2 Systematic uncertainties 

The systematic uncertainties on the peak position and the width, determined from the fit 
shown in Fig. HI were evaluated by changing the selection cuts and the fitting procedure. 
The largest uncertainties on the peak position and the Gaussian width for each check are 
given in parentheses (in MeV). The following systematic studies have been evaluated: 

• the DIS selection cuts were found to have a negligible effect on the peak position. The 
largest uncertainty was found by raising the cut to 30 GeV^ (+0.8, +0.3); 

• the 40-hit requirement for the proton candidates was not used (+0.1, +1.5). The 
variations of the cuts on the tracks in the laboratory frame were found to have a 
negligible effect on the peak; 

• the dE/dx cut was increased to 1.3 mips (+0.6, —0.7). The maximum momentum of 
protons was varied within 1.3 — 4.0 GeV (l^s, -q^)', 

• the bin size was raised and lowered by 1 MeV (loi? +os )! log-likelihood method 
was used for the fitting procedure instead of the method (no change); 

• the fit was made using only a single Gaussian (+0.6, —1.2), or with three Gaussians 
(—0.1, +0.7); the background function was changed to a third-order polynomial (+0.6, 
0.0); 

• the CTD momentum calibration uncertainty on the peak position was calculated from 
the mass measurements of Kg, A and K*, as well as the reconstructed in the Kgp{p) 
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decay channel (IqI)- 

The overall systematic uncertainties of MeV and MeV on the peak position and 
the width were determined by adding the above uncertainties in quadrature. 

9 Summary and conclusions 

The Kg p (p) invariant mass spectrum has been studied in inclusive deep inelastic ep 
scattering for a large range in the photon virtuality. For ^ 10 GeV^ a peak is seen 
around 1520 MeV. 

The peak position, determined from a fit to the mass distribution in the kinematic region 
Q"^ > 20GeV^ is 1521.5 ± 1.5(stat.)+i:7(syst.) MeV, and the measured Gaussian width 
of cr = 6.1 ± 1.6(stat.)^^4(syst.) MeV is above, but consistent with, the experimental 
resolution of 2.0 ± 0.5 MeV. The number of events ascribed to the signal by this fit is 
221 ± 48. The statistical significance, estimated from the number of events assigned to 
the signal by the fit, varies between 3.9a and 4.6(j depending upon the treatment of the 
background. The probability of a similar signal anywhere in the range 1500-1560 MeV 
arising from fluctuations of the background is below 6 x 10~^. 

The results provide further evidence for the existence of a narrow baryon resonance con- 
sistent with the predicted 9+ pentaquark state with a mass close to 1530 MeV and a 
width of less than 15 MeV. Evidence for such a state has been seen by other experiments, 
although the mass reported here lies somewhat below the average mass of these previ- 
ous measurements. In the G"*" interpretation, the signal observed in the Kg p channel 
corresponds to first evidence for an antipentaquark with a quark content of uudds. The 
results, obtained at high energies, constitute first evidence for the production of such a 
state in a kinematic region where hadron production is dominated by fragmentation. 

Acknowledgements 

We thank the DESY Directorate for their strong support and encouragement. The re- 
markable achievements of the HERA machine group were essential for the successful 
completion of this work and are greatly appreciated. We are grateful for the support of 
the DESY computing and network services. The design, construction and installation of 
the ZEUS detector have been made possible owing to the ingenuity and effort of many 
people from DESY and home institutes who are not listed as authors. 



9 



References 

[1] LEPS Coll., T. Nakano et al, Phys. Rev. Lett. 91, 012002 (2003); 
SAPHIR Coll., J. Barth et al., Phys. Lett. B 572, 127 (2003); 
CLAS Coll., V. Kubarovsky et al., Phys. Rev. Lett. 91, 252001 (2003); 
CLAS Coll., V. Kubarovsky et al., Phys. Rev. Lett. 92, 032001 (2004). Erratum; 
ibid, 049902. 

[2] D. Diakonov, V. Petrov and M.V. Polyakov, Z. Phys. A 359, 305 (1997). 

[3] DIANA Coll., V.V. Barmin et al., Phys. Atom. Nucl 66, 1715 (2003); 

A. E. Asratyan, A.G. Dolgolenko, M.A. Kubantsev, Preprint hep-ex/0309042, 2003; 
SVD Coll., A. Aleev et al.. Preprint hep-ex/0401024, 2004; 

HERMES Coll., A. Airapetian et al., Phys. Lett. B 585, 213 (2004); 
COSY-TOF Coll., M. Abdel-Bary et al.. Preprint hep-ex/0403011, 2004. 

[4] NA49 Coll., C. Alt et al., Phys. Rev. Lett. 92, 042003 (2004). 

[5] HI Coll., A. Aktas et al.. Preprint hep-ex/0403017, 2004. 

[6] Particle Data Group, K. Hagiwara et al., Phys. Rev. D 66, 010001 (2002). 

[7] ZEUS Coll., U. Holm (ed.). The ZEUS Detector. Status Report (unpublished), 
DESY (1993), available on http://www-zeus.desy.de/bluebook/bluebook.htinl. 

[8] N. Harnew et al., Nucl. Inst. Meth. A 279, 290 (1989); 

B. Foster et al., Nucl. Phys. Proc. Suppl. B 32, 181 (1993); 
B. Foster et al, Nucl. Inst. Meth. A 338, 254 (1994). 

[9] ZEUS Coll., J. Breitweg et al., Phys. Lett. B 481, 213 (2000); 
ZEUS Coll., J. Breitweg et al., Eur. Phys. J. C 18, 625 (2001). 

[10] M. Derrick et al, Nucl. Inst. Meth. A 309, 77 (1991); 
A. Andresen et al., Nucl. Inst. Meth. A 309, 101 (1991); 
A. Caldwell et al., Nucl. Inst. Meth. A 321, 356 (1992); 
A. Bernstein et al, Nucl. Inst. Meth. A 336, 23 (1993). 

[11] A. Bamberger et al., Nucl. Inst. Meth. A 382, 419 (1996); 

S. Magill and S. Chekanov, Proceedings of the IX Int. Conference on Calorimetry 
(Annecy, Oct 9-14, 2000), B. Aubert et al. (ed.), p. 625. Frascati Physics Series 21, 
Annecy, France (2001). 

[12] A. Bamberger et al, Nucl. Inst. Meth. A 401, 63 (1997). 

[13] J. Andruszkow et al.. Preprint DESY-92-066, DESY, 1992; 
ZEUS Coll., M. Derrick et al., Z. Phys. C 63, 391 (1994); 
J. Andruszkow et al.. Acta Phys. Pol. B 32, 2025 (2001). 



10 



[14] L. Ldnnblad, Comp. Phys. Comm. 71, 15 (1992). 

[15] A. Kwiatkowski, H. Spiesberger and H.-J. Mohring, Comp. Phys. Comm. 

69, 155 (1992). Also in Proc. Workshop Physics at HERA, Ed. W. Buchmiiller and 
G.Ingelman, (DESY, Hamburg, 1991). 

[16] H. Spiesberger, HERACLES anc? D J angoh.- Event Generation for ep Interactions at 
HERA Including Radiative Processes, 1998, available on 
http : //www . desy . de/~hspiesb/dj angoh . html. 

[17] B. Andersson et al., Phys. Rep. 97, 31 (1983). 

[18] T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

[19] R. Brun et al, geantS, Technical Report CERN-DD/EE/84-1, CERN, 1987. 

[20] F. Jacquet and A. Blondel, Proceedings of the Study for an ep Facility for Europe, 
U. Amaldi (ed.), p. 391. Hamburg, Germany (1979). Also in preprint DESY 79/48. 

[21] H. Abramowicz, A. Caldwell and R. Sinkus, Nucl. Inst. Meth. A 365, 508 (1995). 

[22] S. Capstick et al, Phys. Lett. B 570, 185 (2003). 

[23] B. Levchenko, Preprint hep-ph/0401122, 2004. 



11 



Fit 




Gaussian+Bkg. 


2 Gaussians + Bkg. 




M<1700 MeV 


51/47 


35/44 


Peak 1 


mass (MeV) 
width (MeV) 
events 




1465.1 ±2.9 
15.5 ±3.4 
368 ± 121 


Peak 2 


mass (MeV) 
width (MeV) 
events 


1522.2 ± 1.5 
4.9 ± 1.3 

155 ± 40 


1521.5 ± 1.5 
6.1 ± 1.6 

221 ±48 



Table 1: Fit results for > 20 GeV^. 
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Figure 1: The tt+tt" invariant-mass distribution for > \ GeV"^ . The solid 
line shows the fit result using a double Gaussian plus a linear background, while the 
dashed line shows the linear background. 
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Figure 2: The dE/dx distribution as a function of the track momentum for: (a) 
positive and (b) negative tracks. The dE/dx is normalised to a minimum ionising 
particle, defined as the average truncated mean of pion tracks in the momentum 
range 0.3 < p < 0.4 GeF. The solid lines indicate the (anti)proton bands used 
in this analysis, and the dotted line denotes the cuts dE/dx > 1.15 mips and 
p < 1.5 GeV. 
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Figure 3: Invariant-mass spectrum for the Kg p (p) channel, after the cuts f < 
dE/dx < F, dE/dx > 1.15 mips andp < 1.5 GeV , integrated for (a) Q"^ > 1 GeV"^ , 
(h) Q2 > 10 GeV^, (c) g2 > 30 GeV^ , (d) > 50 GeV\ (e) > 1 GeV^ and 
W < 125 GeV (f) > i QeV"^ and W > 125 GeV. The histogram shows the 
prediction of the Ariadne MG simulation normalised to the data in the mass 
region above 1650 MeF. 
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Figure 4: Invariant-mass spectrum for the K^p (p) channel for > 20 GeV'^ , 
with other cuts as in Fig. The solid line is the result of a fit to the data using 
a three-parameter background function plus two Gaussians (see text). The dashed 
lines show the Gaussian components and the dotted line the background according 
to this fit. The histogram shows the prediction of the Ariadne MG simulation 
normalised to the data in the mass region above 1650 MeF. The inset shows the 
Kgp (open circles) and the Kg p (black dots) candidates separately, compared to 
the result of the fit to the combined sample scaled by a factor of 0. 5. 
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